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We study magnetic field-induced three-dimensional ordering transitions in low-dimensional quan- 
tum spin liquids, such as weakly coupled, antiferromagnetic spin-1/2 Heisenberg dimers and ladders. 
Using stochastic series expansion quantum Monte Carlo simulations, thermodynamic response func- 
tions are obtained down to ultra-low temperatures. We extract the critical scaling exponents which 
dictate the power-law dependence of the transition temperature on the applied magnetic field. These 
are compared with recent experiments on candidate materials and with predictions for the Bose- 
Einstein condensation of magnons obtained in mean-field theory. 



PACS numbers; 75.10.Jm, 75.30.-m, 75.50.Ee 

Many low-dimensional quantum spin liquids, such as 
antiferromagnetic Heisenberg spin ladders, have a non- 
magnetic valence-bond ground state with a finite energy 
gap to the lowest band of triplet excitations. This spin 
gap can be suppressed by an applied magnetic field. In- 
creasing the external field beyond a critical value hd 
leads to partial spin-polarization and incommensurate, 
gapless excitations. A second transition occurs at a 
higher critical field, hc2, above which the system becomes 
fully polarized. It has been suggested that in the par- 
tially polarized phase, hci < h < hc2, weak couplings 
between the low-dimensional subsystems induce a three- 
dimensional (3D) condensation of magnetic excitations 
at low temperatures. 

Indications of such 3D ordering transitions induced by 
external magnetic fields were observed by high-field nu- 
clear magnetic resonance and inelastic neutron scatter- 
ing measurements on the quantum spin liquids TlCuCla 
§1 and Cu2(C5Hi2N2)2Cl4 The crystal struc- 

ture of TlCuCla suggests that this compound consists 
of weakly coupled CU2CI6 dimers, involving two neigh- 
boring spin- 1/2 Cu2+ ions in its b-c plane. |jl|] A 
small spin gap of A « 7.5K was determined from mea- 
surements of the magnetic susceptibility. [p|JTT| On the 
other hand, Cu2(C5Hi2N2)2Cl4 is generally viewed as a 
compound of weakly coupled two-leg ladders with intra- 
ladder exchange constants J± w 13. 2K, J|| sa 2.5K, 
and a spin gap A « 10.5K. However, a very recent 
study has concluded that this material consists of frus- 
trated planar networks of dimers. We may further 
expect field-induced 3D ordering phenomena to occur in 
(C5Hi2N)2CuBr4. This compound is believed to con- 
sists of weakly coupled two-leg ladders with a strong 
rung coupling, J± ~ 13. 3K, between adjacent CuBr4 
tetrahedra, and a weaker coupling, J|| ~ 3.8K along the 
legs of the ladders. |0 Its spin gap was observed to be 
A w 9.5K. 

Whatever the effective dimensionalities of these com- 
pounds may turn out to be (d=0 for TlCuCls, 
d=l or d=2 for Cu2(C5Hi2N2)2Cl4, and d^l for 



(C5Hi2N)2CuBr4), they share the essential common fea- 
ture of a spin-liquid ground state and a small but finite 
spin gap that can be overcome by presently accessible 
magnetic fields. 

From the theoretical side, it has been proposed that the 
observed field-induced 3D magnetic ordering transition 
in these systems can be interpreted as a Bose-Einstein 
condensation of magnons, resulting in a staggered, trans- 
verse magnetic order. ||l|,|l^-jl^ 

This class of quantum phase transition is character- 
ized by a critical exponent a, which relates the ordering 
temperature Tc to the applied magnetic field according 
to 
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The exact determination of this power-law dependence 
from experimental data turns out to be somewhat deli- 
cate, because it requires simultaneous fitting of the crit- 
ical field and the scaling exponent. Furthermore, the 
number of data points tends to be sparse in the scaling 
regime around {hc,Tc). Nevertheless, the experimentally 
deduced exponents, a « 2.0 for TlCuCla and a ~ 1.5 
for Cu2(C5Hi2N2)2Cl4, appear to be within the general 
range of the value ajBEC — 3/2, expected from isotropic 
Bose-Einstein Hartree-Fock theory. ||,|l^,|l^ However, 
material-specific details such as spin-phonon coupling 
and magnetic frustration in these compounds are likely 
to influence the precise values of the scaling exponents. 

In this paper, we address the fundamental question of 
obtaining these critical scaling properties directly from 
microscopic models of weakly coupled low-dimensional 
quantum spin systems. We apply the recently devel- 
oped stochastic series expansion quantum Monte Carlo 
(QMC) technique to the 3D antiferromagnetic spin- 
1/2 Heisenberg model with spatially anisotropic ex- 
change couplings. 
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where h denotes the applied magnetic field. The rela- 
tive strengths of the nearest-neighbor exchange coupling 
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constants Jij are illustrated in Fig. 1, in which planar sec- 
tions of two clusters with different configurations of Jij 
are shown. Fig. 1(a) illustrates an ensemble of weakly 
coupled dimers oriented along the x-direction, a con- 
figuration resembling the minimal effective magnetically 
structure of TlCuCla. Fig. 1(b) shows a quasi-lD array of 
weakly coupled two-leg Heisenberg ladders oriented along 
the y-direction, as may be realized in (C5Hi2N)2CuBr4 

or Cu2(C5Hi2N2)2Cl4. 
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FIG. 1. Layers of anisotropically coupled quantum spins: 
(a) weakly coupled dimers, and (b) weakly coupled ladders. 
The strength of the exchange coupling constants is indicated 
by the thickness of the lines. The 3D inter-layer couplings (not 
shown in the figure) have the same strength as the weakest 
links (dashed lines) within the planes. 

The numerical algorithm involves expansions of the 
partition function in inverse temperature, uses local and 
global system updates, and is significantly more efficient 
than conventional QMC schemes. In this work we study 
cubic lattices with up to 10 x 10 x 10 quantum spins down 
to very low temperatures, by which is meant less than 1% 
of the exchange coupling strength J which sets the scale 
of the problem. Moreover, the stochastic series expan- 
sion QMC method can handle external magnetic fields of 
any strength without the problems common to worldline 
QMC techniques, such as the QMC loop algorithm. 

In Fig. 2(a), the low-temperature regime of the uniform 
magnetization is shown for the system of weakly coupled 
dimers depicted in Fig. 1(a) at various magnetic fields 
within the partially polarized regime, hd < h < hc2- The 
anisotropy of the exchange coupling constants, J' / J — 
1/15, was chosen according to estimates for TlCuCla. 
Here J is the strong intra-dimer and J' the weak inter- 
dimer coupling. In this case, the observed critical fields 
are well approximated by perturbation theory about the 
limit of non- interacting dimers, giving hd — J — 5 J'/2 
and hc2 — J + 5 J'. At ultra-low temperatures on the or- 
der of J', the magnetization curves of the weakly coupled 
dimers (solid lines) have a maximum if h„i < h < hc2 and 
a minimum if hd < h < /i„j, where hm = {hd + /ic2)/2, 
indicating the onset of a magnetic field-induced 3D order- 
ing. This feature, emphasized by large filled circles, is ab- 
sent in the magnetization curves for the non-interacting 



limit (J' = 0), denoted by the dashed lines. While effects 
of the weak inter-dimer couplings are clearly negligible 
at higher temperatures, their relevance at low temper- 
atures is seen in the departure of the solid lines from 
the dashed lines. In contrast to the mean-field theory 
of weakly coupled dimers, p9[ the magnetization is ob- 
served to be strongly temperature-dependent below the 
ordering temperature Tc, as shown in Figs. 2(b) and (c). 
On the other hand, the 3D behavior is in good qualitative 
agreement with the recently proposed Bose-Einstein con- 
densation (EEC) description of this ordering transition. 




0.42 


(b) 




0.41 


— ' — 1 — ' — 1 








' ^ 


0.09 


(c) 




0.08 







0.0 0.1 0.2 

T/J 



0.3 0.00 0.05 0.10 
T/J 



FIG. 2. (a) Temperature dependence of the uniform mag- 
netization in the system of weakly coupled dimers shown in 
Fig. 1(a). The coupling anisotropy is J' /J = 1/15, and the 
magnetic fields are h/J = 0.80, 0.86, 0.90, 0.97, 1.08, 1.17, 
1.23, 1.27, 1.32, 1.37. Extrema in m{T) indicate the onset 
of 3D ordering, and are denoted by enlarged filled circles. 
Results for the non-interacting {J' /J = 0) limit are plotted 
with dashed lines, (b) Low-temperature regime of m{T) at 
h = 1.27J and (c) at h = 0.9J. 

The dependence of Tc on the applied magnetic field can 
be extracted from the locus of the extrema in m{T) which 
are robust to finite size effects for the system sizes consid- 
ered here. The resulting phase diagram is plotted in 
Fig. 3 for the ensemble of weakly coupled dimers. As ex- 
pected from Eq. 1, the transition temperature exhibits a 
power-law dependence in the vicinity of hd and hc2 ■ Our 
best fits, shown in the insets of Fig. 3, yield a = 2.7±0.2 
for the lower-field transition and a = 2.3 ± 0.2 in the 
vicinity of hc2- These critical exponents differ signifi- 
cantly from the value obtained by standard, isotropic 
Bose-Einstein Hartree-Fock theory, ajBEC = 3/2. It 
was previously pointed out that some of this deviation 
may be attributed to shortcomings of the Hartree-Fock 
description in the critical regions. |15 



Another important point is that the dispersion of the 
spin triplet band is strongly anisotropic. pi| For the case 



2 



of linearly aligned, weakly coupled dimers the dominant 
feature in the triplet excitation spectrum is a parabolic 
bonding band along the strong-coupling direction which 
becomes populated when the magnetic field is raised be- 
yond hci- Q Along the weak-coupling directions the 
triplet excitation spectrum is linear in this regime. Fol- 
lowing the BEC Hartree-Fock treatment of , a crit- 
ical exponent a = 5/2 is obtained for this case. ||2^] 
This anisotropic BEC Hartree-Fock exponent is in better 
agreement with the QMC simulations than ajsEC = 3/2. 
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second feature, marked by the filled arrows in Figs. 4(b) 
and (d), occurs at still lower temperature. This is the 3D 
magnetic ordering transition, also manifested as a change 
of slope in the energy, as shown in the insets. Measure- 
ments of the magnetic specific heat in Cu2(C5Hi2N2)2Cl4 
indeed show such 3D ordering features at temperatures 
below the onset of quantum critical behavior. |q| 
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FIG. 3. Phase diagram of the magnetic field-induced or- 
dering transition in a system of weakly coupled dimers. The 
scaling exponents in the vicinity of the lower and upper criti- 
cal fields are extracted from double- logarithmic plots of Tc{h), 
shown in the insets. 

On the other hand, measurements of the lower-field 
critical exponent in TlCuCla give a ~ 2.0 H) and a ~ 2.2 
^ . This suggests that the generic Heisenberg Hamilto- 
nian which we have studied is not sufficient to accurately 
model the critical properties of this compound. Reasons 
for this discrepancy may be found in (i) the much more 
complex bandstructure of the spin-triplet excitations in 
the real material, |^ , ^ (ii) magnetic frustration effects, 
and (iii) coupling to phonon degrees of freedom. 

Let us now turn to the system of weakly coupled two- 
leg ladders depicted in Fig. 1(b). For the couphng con- 
stants chosen in Fig. 4 and 5, the critical magnetic fields 
are hd = 0.61J and hc2 — 1.86J. It has recently been 
pointed out that this system exhibits extrema in m(T) 
even in the complete absence of inter-ladder couplings. 
p5[ We can confirm this observation with our simulations 
of the uniform magnetization, as shown in Figs. 4(a) and 
(c) . At magnetic fields slightly above hci minima occur at 
low temperatures, whereas maxima are observed at larger 
magnetic fields close to /ic2- These features indicate the 
transition into low-temperature Luttinger Liquid behav- 
ior in the partially polarized regime, hd < h < hc2- |2^] 
When the weak inter-ladder couplings are switched on, a 



FIG. 4. Low-temperature regime of the uniform magneti- 
zation in ladder systems, (a) and (c): isolated two- leg lad- 
ders, (b) and (d): system of weakly coupled ladders shown 
in Fig. 1(b). The coupling constants are J = J_l, Jy — J/3, 
J' = for (a) and (c), and J' = J/15 for (b) and (d). The 
applied magnetic fields are h = 0.83 J for (a) and (b), and 
h — 1.5J for (c) and (d). The temperature dependence of the 
energy is shown in the insets. 

A rich magnetic phase diagram, shown in Fig. 5, can be 
constructed from these magnetization response functions, 
and shares some of the essential features reported by ex- 
periments on (C5Hi2N)2CuBr4 and Cu2(C5Hi2N2)2Cl4. 
At high temperatures (not shown here), the system is 
effectively zero-dimensional, with a paramagnetic Curie- 
law temperature dependence of the magnetic susceptibil- 
ity. Lowering the temperature, there is a low-field disor- 
dered spin liquid regime {h < h„i) and a high-field spin- 
polarized phase {h > hm), where hm ~ {hd + /ic2)/2 
separates these two regions. The onset of Luttinger 
Liquid behavior is found in the partially spin-polarized 
regime, hd < h < hc2, at temperatures below 0.13 J for 
J|| = J/3. The effective dimensionality of this region is 
d=l. Finally, 3D ordering occurs at still lower temper- 
atures, as indicated by the lowest transition line in Fig. 
5. The scaling exponents extracted from fits of Tc{hc) 
to the power law in Eq. 1 are a = 3.1 ± 0.2 at the 
lower critical field and a = 1.8 ± 0.2 at the upper criti- 
cal field. These values are quite different from isotropic 
Bose-Einstein Hartree-Fock theory, suggesting that a suc- 
cessful effective theory needs to account more accurately 
for the quantum dynamics of the lower-dimensional sub- 
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FIG. 5. Magnetic phase diagram of weakly coupled ladders. 
The upper lines indicate the transition into a Luttinger Liquid 
regime, also present in isolated two- leg ladders. The lower 
line marks the magnetic field induced 3D ordering transition. 
The scaling exponents in the vicinity of the lower and upper 
critical fields are extracted from double-logarithmic plots of 
Tc{h), shown in the insets. 

In summary, we have studied the magnetic phase dia- 
gram and the critical behavior in models of weakly cou- 
pled quantum spin liquids. Using the stochastic series 
expansion QMC method, we were able to reach suffi- 
ciently low temperatures that magnetic field-induced 3D 
ordering could be observed. The scaling exponents de- 
pend strongly on the dimensionality and on the quantum 
dynamics of the subsystems, Heisenberg dimers and lad- 
ders, which reflect the strongly anisotropic dispersion of 
the triplet excitation bands. While this study has con- 
centrated on the fundamental features of generic model 
Hamiltonians, much more detailed models are clearly 
needed to account for the specific scaling properties ob- 
served in real materials. 
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